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It is shown that photon shot noise and radiation-pressure back-action noise are the sole forms 
of quantum noise in interferometric gravitational wave detectors that operate near or below the 
standard quantum limit, if one filters the interferometer output appropriately. No additional noise 
arises from the test masses' initial quantum state or from reduction of the test-mass state due to 
measurement of the interferometer output or from the uncertainty principle associated with the test- 
mass state. Two features of interferometers are central to these conclusions: (i) The interferometer 
output (the photon number flux TV (t) entering the flnal photodetector) commutes with itself at 
different times in the Heisenberg Picture, [M{t),JV{t')] = and thus can be regarded as classical, 
(ii) This number flux is linear to high accuracy in the test-mass initial position and momentum 
operators Xo and po, and those operators influence the measured photon flux Jy{t) in manners that 
can easily be removed by filtering. For example, in most interferometers Xo and po appear in Af{t) 
only at the test masses' ~ 1 Hz pendular swinging freqency and their infiuence is removed when the 
output data are high-pass filtered to get rid of noise below ~ 10 Hz. 

The test-mass operators Xo and Po contained in the unfiltered output Af{t) make a nonzero 
contribution to the commutator [M{t),J\f{t')]. That contribution is precisely cancelled by a nonzero 
commutation of the photon shot noise and radiation-pressure noise, which also are contained in 
J\(t). This cancellation of commutators is responsible for the fact that it is possible to derive 
an interferometer's standard quantum limit from test-mass considerations, and independently from 
photon-noise considerations, and get identically the same result. 

These conclusions are all true for a far wider class of measurements than just gravitational-wave 
interferometers. To elucidate them, this paper presents a series of idealized thought experiments 
that are free from the complexities of real measuring systems. 
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I. QUESTIONS TO BE ANALYZED AND 
SUMMARY OF ANSWERS 

It has long been known that the Hciscnbcrg uncer- 
tainty principle imposes a "standard quantum limit" 
(SQL) on high-precision measurements [1-3]. This SQL 
can be circumvented by using "quantum nondemolition" 
(QND) techniques [2 9]. 

For broad-band interferomctric gravitational-wave de- 
tectors the SQL is a limiting (single-sided) spectral den- 
sity 



= m/4 



(1.3) 



Shif) 



m(27r/)2L2 



(1.1) 



for the gravitational- wave field h{t) [10,11]. Here fi is 
Planck's constant divided by 2it, m, is the mass of each 
of the interferometer's four test masses, L is the interfer- 
ometer's arm length, and / is frequency. 

This SQL firmly constrains the sensitivity of all con- 
ventional interferometers (interferometers with the same 
optical topology as LIGO's first-generation gravitational- 
wave detectors) [12,13]. LIGO's second-generation inter- 
ferometers (LIGO-II; ca. 2008) are expected to reach this 
SQL for their m = 40 kg test masses in the vicinity of 
/ ~ 100 Hz [14], and may even beat it by a modest 
amount thanks to a "signal recycling mirror" that con- 
verts them from conventional interferometers into QND 
devices [15-17]. LIGO-III interferometers are likely to 
beat the SQL by a factor ~ 4 or more; see, e.g., [13]. 

In the research and development for LIGO-II interfer- 
ometers [14-17] and in the attempts to invent strongly 
QND LIGO-III interferometers [18-24,13], it is important 
to understand clearly the physical nature of the quantum 
noise which imposes the SQL, and to be able to com- 
pute with confidence the spectral density of this quantum 
noise for various interferometer designs. These issues are 
the subject of this paper. 

There are two standard ways to derive the 
gravitational- wave SQL (1.1), and correspondingly two 
diff'erent viewpoints on it. The first derivation [10,25] fo- 
cuses on the quantum mechanics of the interferometer's 
test masses and ignores the interferometer's other details. 
In the simplest version of this derivation, one imagines a 
sequence of instantaneous measurements of the difference 



x= {xi- X2) - {xs - X4) 



(1.2) 



of the center-of-mass positions of the four test masses, 
and from this measurement sequence one infers the 
changes of x and thence the time varying gravitational- 
wave field h{t) = x{t)/L. At time t immediately after 
one of the measurements, the test masses' reduced state 
has position variance [Ax(t)]2 no smaller than the mea- 
surement's accuracy. During the time interval t = t' — t 
between this measurement and the next, the test masses 
are free, so x{t) evolves as the position of a free particle 
with mass 



[the reduced mass of the four-body system with relative 
position (1.2)]. The Heisenberg-Picture commutation re- 
lations for a free particle 



ih{t' - t) AihT 



m 



(1.4) 



imply that, whatever may be the state of the test masses, 

the variance [/\x[t')Y of x just before the next measure- 
ment must satisfy the Heisenberg uncertainty relation 



\x{t)\x{t') > 



n\t - 1'\ _ 2hT 

2/i m 



(1.5) 



The accuracy with which the change of x between t 
and t' can be measured is no better than the value ob- 
tained by setting Ax{t) = A.T(i'), and in classical lan- 
guage that accuracy is related to the minimum possi- 
ble spectral density of the noise at frequency / ~ l/nr 
by Ax{t) = Ax{t') ~ \/Sh{f)/T. Simple algebra then 
gives expression (1.1) for the SQL of Sh{f)- A more 
sophisticated analysis [10], based on measurements that 
are continuous rather than discrete and on a nonunitary 
Feynman-patli-integral evolution of the test-mass state 
[26,27], gives precisely the SQL (1.1). 

The second derivation of the SQL [28,29] ignores the 
quantum mechanics of the test mass, and focuses instead 
on that of the laser light which monitors the test-mass 
motion. The light produces two kinds of noise: pho- 
ton shot noise, which gets superposed on the output 
gravitational-wave signal, and radiation-pressure fluctu- 
ations, which produce a random back-action force on the 
test masses, thereby influencing their position evolution 
and thence the interferometer output. In an ideal, SQL- 
limited interferometer, both noises — shot and radiation- 
pressure — arise from quantum electrodynamic vacuum 
fluctuations that enter the interferometer through its 
dark port and superpose on the highly classical laser light 
[28,29]. The radiation-pressure spectral density is pro- 
portional to the laser-light power P, the shot-noise spec- 
tral density is proportional to l/P, and their product is 
independent of P and is constrained by the uncertainty 
principle for light (or equivalently by the electromagnetic 
field commutation relations) to be no smaller than 



(1.6) 



[cf. Eqs. (6.7) and (6.17) of [3] in which there is a factor 
1 /4 on the right side because Ref. [3] uses a double-sided 
spectral density, while the present paper uses the gravity- 
wave community's single-sided convention]. In Eq. (1.6) 
Sxif) is the spectral density of the shot noise that is 
superposed on the interferometer's output position sig- 
nal x{t), Spif) is the spectral density of the radiation- 
pressure force that acts on the test-mass ccnter-of-mass 
degree of freedom x, and we have assumed that the shot 
noise and radiation-pressure force are uncorrected as is 
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the case for conventional (LIGO-I type) interferometers 
[13,15-17]. At frequency / the test mass responds to 
the Fourier component F{f) of the force with a position 
change x{f) = — F(/)/[/i(27r/)^], and correspondingly 
the net gravitational-wave noise is 

W) = ^(^. + ^^) . (1.7) 

By combining Eqs. (1.6), (1.7) and (1.3), we obtain the 
SQL (1.1) for a conventional interferometer, e.g. LIGO-L 
In view of these two very different derivations of 
the SQL, test-mass quantization and light quantization, 
three quc!stions arise: (i) Arc the test-mass quantization 
and the light quantization just two different viewpoints 
on the same physics? — in which case the correct SQL 
is Eq. (1.1). Or arc they fully or partially independent 
effects? — in which case we would expect their noises to 
add, causing the true SQL for Sh to be larger by, perhaps, 
a factor 2 (and thence the event rate in an SQL-limited 
interferometer to be reduced by a factor ~ (V^)"^ — 3). 

(ii) How should one compute the quantum noise in can- 
didate designs for the QND LIGO-II and LIGO-III inter- 
ferometers? One inevitably must pay close attention to 
the behavior of the light (and thus also its quantization), 
since the optical configuration will differ markedly from 
one candidate design to another. Must one also pay close 
attention to the quantum mechanics of the test masses, 
including their commutation relation (1.4) and the con- 
tinual reduction of their state as information about them 
is continually put onto the light's modulations and then 
measured? (iii) Similarly, how should one design a QND 
interferometer? Need one adjust one's design so as to 
drive both the light's noise and the test-mass noise be- 
low the SQL? 

As we shall show, the answers are these: (ii) The test- 
mass quantization is irrelevant to the interferometer's 
noise and correspondingly test-mass state reduction is 
irrelevant, if one filters the output data appropriately. 
(For interferometers with conventional optical topology 
such as LIGO-I, it is sufficient to discard all data near 
the test masses' ^ 1 Hz swinging frequency.) Therefore, 
one can ignore test-mass quantization and state reduction 
when computing the noise of a candidate interferometer. 

(iii) Similarly, one can ignore the test mass's quantum 
noise when designing a QND interferometer that beats 
the SQL. One need only pay attention to the light's quan- 
tum noise, and in principle, by manipulating the light ap- 
propriately (and filtering the output data appropriately), 
one can circumvent the SQL completely, (i) Correspond- 
ingly, the SQL (1.1) as derived from light quantization 
is precisely correct; there is no extra factor 2 caused by 
test-mass quantization. [The fact that one can also derive 
the SQL from test-mass quantization is a result of an in- 
timate connection between the uncertainty principles for 
a measured system (the test masses in our case) and the 
system that makes the measurement (the light). We shall 
elucidate this intimate connection from one viewpoint at 



the end of Sec. TIB 4. From another viewpoint, it is due 
to the fact that the commutator [x(t),x{t')], which un- 
derlies the test-mass derivation (1.4), (1.5) of the SQL, 
also underlies the derivation of the measuring light's un- 
certainty relation (1.6); see the role of the generalized 
susceptibility x(*)*') = {^/i^)[x{t'),x{t)] in Sec. 6.3 of 
Ref. [3].] 

Central to our answers (i), (ii) and (iii) is the fact 
that an interferometric gravitational-wave detector does 
not monitor the time-evolving test-mass position x{t). 
Rather, it only monitors classical changes in x{t) in- 
duced by the classical gravitational- wave field h{t) and 
other classical^ forces (thermal, seismic, ...) acting on 
the test masses, and it does so without extracting infor- 
mation about the actual quantized position x{t). The 
detector has a classical input [h{t)] and a classical out- 
put [h{t) contaminated by noise that (as we shall see) 
commutes with itself at different times and that there- 
fore can be regarded as a time-evolving c- number]. The 
quantum properties of the test masses and the light are 
merely intermediaries through which the classical signal 
must pass. This would not be the case for a device de- 
signed to make a sequence of absolute measurements of 
the quantum mechanical position x{t). 

Our answers (i), (ii), (iii) hold true for a far wider 
range of measuring devices than just interferometric 
gravitational- wave detectors. They hold quite generally 
for any well-designed device that measures a classical 
force acting on any quantum mechanical system. In par- 
ticular, they remain true if the device makes measure- 
ments that are linear in the sense of Appendix B, and 
one filters the device's output to remove all information 
at the natural frequencies of the quantum system's dy- 
namics (e.g. at its eigenfrequency if the quantum system 
is a harmonic oscillator). 

In Sec. II we will elucidate these answers by consider- 
ing pedagogical examples of idealized devices that make 
discrete, quick measurements on a test mass. These ex- 
amples will reveal two central underpinnings of our an- 
swers: (a) the vanishing of the measurement's "output 
commutators" — i.e., the commutators of the observ- 
ables (Hermitian operators) that represent the entries in 
the output data stream, and (b) a data-processing pro- 
cedure that removes from the data all influence of the 
test-mass quantum observables (initial position Xq and 
initial momentum po). Our examples will also elucidate 
two strategies for beating the SQL: (A) put the mea- 
suring apparatus ("meters") into specially chosen initial 
states (the analog of squeezed states), and (B) measure 
a wisely chosen linear combination of position and mo- 



^All these forces — gravitational- wave, thermal, seismic, etc. 
— actually do have a quantum component, but in practice 
their levels of excitation are so large that we can regard them 
as classical. 
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mcntiini for the test mass and thereby remove the effects 
of the meters' back action from the output data (make a 
"quantum variational measurement" ) . 

Our examples are the following: We will begin in Sec. 
II A with a simple, idealized, instantaneous single mea- 
surement of the position of a single test mass. This ex- 
ample will demonstrate that the noise associated with 
test-mass quantization and the noise associated with 
the meter's quantization are truly independent (though 
closely linked), and will illustrate how under some cir- 
cumstances they can add, producing a doubling of the 
noise power. Then, in Sec. TIB, we will analyze the use 
of a sequence of these idealized, instantaneous position 
measurements to monitor a classical force that acts on 
the test mass. This example will illustrate the vanishing 
self-commutator of the output data samples, which arises 
from a cancellation of the test-mass-position commuta- 
tor by the measurement-noise commutator; it also will 
illustrate how signal processing can remove all influence 
of test-mass quantization and test-mass state reduction 
from the output data stream. Our third example (Sec. 
II C) will be a Heisenberg-microscope-like realization of 
these instantaneous, idealized position measurements, in 
which a pulse of near- monochromatic light is reflected off 
the test mass, thereby encoding the test-mass position in 
a phase shift of the light. This example will give real- 
ity to the idealized examples in Sees. II A and II B, and 
will help connect them to the subsequent discussion of 
interferometric gravitational-wave detectors. 

In Sec. Ill we will use the insights from our peda- 
gogical examples to prove and elucidate our three an- 
swers [(i), (ii), (iii) above] for gravitational- wave inter- 
ferometers, and also for a wide range of other classical 
force measiircmcnts. The underpinnings for our answers 
will be: (a) a proof that for a quantized electromagnetic 
wave, such as that entering the final photodetector of an 
interferometer, the photon number flux operator com- 
mutes with itself at different times (this flux is the out- 
put data stream), and (b) a proof that all influence of 
the test-mass quantum observables can be removed from 
the output data stream by appropriate flltering, and for 
conventional interferometers it is sufficient to remove all 
data near the test masses' ~ 1 Hz swinging frequency, 
e.g. by the kind of high-pass filtering that is routinely 
used in gravitational-wave detectors. Our analysis will 
also elucidate QND interferometer designs based on (A) 
squeezed-input states for light and (B) variational-output 
measurements. 

The issues studied in this paper are most efficiently an- 
alyzed in the Heisenberg picture, and the Hciscnbcrg pic- 
ture gives particularly clear insights into them. For this 
reason, we will use the Heisenberg picture throughout the 
body of this paper. Readers who arc uncomfortable with 
the Heisenberg picture may find Appendix A reassuring; 
there we will give a detailed Schroedinger-picture analy- 
sis of the most important of our pedagogical examples, 
that of Sec. II B 



II. PEDAGOGICAL EXAMPLES 

A. A single position measurement: 
"double" uncertainty relation 

We begin with a simple pedagogical example of a single 

measurement of the position of a single test mass. The 
Heisenberg microscope is a famous realization of this ex- 
ample; see Sec. II C. 

The measurement is idealized as instantaneous and as 
occuring at time t — i). At times arbitrarily close tot ~ 0, 
the Hamiltonian for the test mass (with position and mo- 
mentum X and p) and the mcasiiring device (the m,eter, 
with generalized position Q and generalized momentum 
P) is 



H = 



El 

2fj, 



5{t)xP - 



2M 



(2.1) 



Here 5{t) is the Dirac delta function, /i is the test mass's 
mass and M is the generalized mass of the meter. For 
pedagogical simplicity we make M arbitrarily large so Q 
and P do not evolve in the Heisenberg Picture except 
at the moment of interaction, and correspondingly we 
rewrite the Hamiltonian as 



H=E-- 5(t)xP . 
2/i ^ ^ 



(2.2) 



A simple calculation in the Heisenberg picture gives 
the following expressions for the positions and momenta 
immediately after the measurement, in terms of those 
immediately before: 



-^after — -^before 5 

»^after — ^before ■> 

Q&iter — Qbcforc '^before 

Paftcr — J^bd'orc + Pb cfore ■ 



(2.3a) 
(2.3b) 
(2.3c) 
(2.3d) 



The meter's generalized position Qafter is amplified and 
read out classically immediately after the interaction, to 
determine the test-mass position. The resulting mea- 
sured position, expressed as an operator, is Xmeas = 
-Qafter = ^before " Qbefore [Eq. (2.3c)], and the mea- 
surement leaves the actual test-mass position operator 
unperturbed [Eq. (2.3b)] but it perturbs the test-mass 
momentum [Eq. (2.3d)]. 

It is instructive to rewrite Eqs. (2.3c) and (2.3d) in the 
form 



•^meas — •^before ~t~ ^•^meas ? 
Pafter = Pbefore + ^PBA , 



with 



Sx„ 



-Q 



bcfo 



Sp 



BA 



(2.4a) 
(2.4b) 



(2.5) 



The simple equations (2.4) embody the measurement re- 
sult and its back action; ^meas is the measured value of 
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^before = iaftcr; (^^mcas IS thc noisc superposcd On that 
measured value by the meter, and SpBA is the back-action 
impulse given to the test mass by thc meter. Equations 
(2.4) are actually much more general than our simple 
example; they apply to any sufficiently quick,^ "linear" 
measurement; see Eqs. (5.2), (5.14) and (5.23) of Ref. [3], 
and see Appendix B below. 

The initial test-mass position and momentum and the 
initial meter position and momentum have the usual com- 
mutation relations, 

[*^before ? Pbefore] — — [^before ^ -^before] • (^*^) 

The second of these and Eqs. (2.5) imply that thc mea- 
surement noise Sxmeas and the back-action impulse 6pBA 
have this same standard commutator, but with the sign 
reversed 

[J^meas, Spba] = -i^ ■ (2.7) 

This has an important implication: The measured value 
of thc test-mass position and the final value of the test- 
mass momentum commute: 

[^meas,Pafter] = . (2.8) 

This result, like thc simple measurement and back-action 
equations (2.4), is true not only for this pedagogical ex- 
ample, but also for any other sufficiently quick, linear 
measurement; see, e.g., Sec. II C below. 

It is evident from Eqs. (2.4) and (2.5) that the vari- 
ances of imcas and Pafter are influenced by the initial 
states of both the meter and the test mass: 

(Aa; 

me as 

f = (AXbefore)^ + (AQbefore)^ , (2.9) 

(Apafter)^ = (Apbefore)^ + (APbefore)^ • (2.10) 

Here wc have assumed, as is easy to arrange, that the 
initial states of the meter and the test mass are uncorre- 
lated. Now, the initial states of the test mass and meter 
are constrained by the uncertainty relations 

AXbefore • Apbefore > ^ , (2-11) 
AQbefore ' APbefore > 2 ' (2-12) 

which follow from thc commutators (2.6). From the 
viewpoint of the measurement equations (2.4), the me- 
ter equation (2.12) is an uncertainty relation between 
the noise (J^moas = —Qbcforc that the meter superim- 
poses on the output signal, and the back-action impulse 
^Pba = -Pbefore that the meter gives to the test mass. In 



^i.c, quick compared to tlic evolution of thc wave function 
of the measured quantity, so it can be regarded as constant 
during the measurement. 



the Heisenberg microscope, Sx^^^as would be photon shot 
noise and SpBA would be radiation-pressure impulse. 

The test-mass uncertainty relation (2.11) and meter 
uncertainty relation (2.12) both constrain the product 
of the measurement error (2.9) and the final momentum 
uncertainty (2.10), and by equal amounts. The result is 
a "doubling" of the uncertainty relation, so 

Aa^meas • Apafter > ' ^ ■ (2.13) 

This doubling of the uncertainty relation relies cru- 
cially on our assumption that the initial states of the 
test mass and meter arc uncorrclatcd. Correlations can 
produce a violation of thc uncertainty relation (2.13). 
For example, initial correlations can be arranged so as 
to produce (in principle) a vanishing total measurement 
error A.Tnicas = and a finite Apafter so the product 
A.Tnicas ' Apaitcr vanishcs — a result permitted by the 
vanishing commutator (2.8). 

B. Monitoring a classical force: 
"single" uncertainty relation 

As we emphasized in Sec. I, the goal of LIGO-type de- 
tectors is not to measure any observables of a test mass, 
but rather to monitor an external force that acts on it. 
Correspondingly, it is desirable to design the measure- 
ment so the output is devoid of any information about 
the test mass's initial state. As we shall see, this is readily 
done in a way that removes the initial-state information 
during data processing. Thc result is a "single" uncer- 
tainty relation: the measurement result is influenced only 
by the quantum properties of the meter and not by those 
of the test mass. 



1. Von Neumann's thought experiment 

We illustrate this by a variant of a thought experiment 

devised by von Neumann [30] and often used to illustrate 
issues in thc quantum theory of measurement; see, e.g., 
[31] and references therein. We analyze this thought ex- 
periment using the Heisenberg picture in the body of 
this paper, and we give a Schroedinger-picture analysis 
in Appendix A. 

Our von Neumann thought experiment is a simple gen- 
eralization of the position measurement described above. 
Specifically, wc consider a free test mass, with mass /i, 
position X and momentum p, on which acts a classical 
force F(t). To monitor F{t), we probe the test mass 
instantaneously at times t = 0, t, .... (N — l)r using 
N independent meters labeled r = 0,1,...,A^ — 1. Each 
meter is prepared in a carefully chosen state, it then in- 
teracts with the test mass, and then is measured. We 
filter the measurement results to deduce F{t). Meter r 
has generalized coordinate and momentum Qr and Pr, 
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and its free Hamiltonian is vanishingly small, so Qr and 
Pr do not evolve except at the moment of interaction. 
The total Hamilton for test mass plus classical force plus 
meters is 



H 



2/x 



N-l 

E 

r=0 



5{t — rT)xPr 



(2.14) 



We denote by xq and po the test-mass position and 
momentum at time t = Q when the experiment begins, 
and by Xr and Pr their values immediately after inter- 
acting with meter r, at time t = rr. The momentum of 
meter r is a constant of the motion, so we denote it by 
Pr at all times. The meter coordinate changes due to the 
interaction; we denote its value before the interaction by 
Qboforc .^-^^ after the interaction by Qr- 

It is easy to show, from the Heisenberg equations for 
the Hamiltonian (2.14), that the test-mass position im- 
mediately after its r'th interaction is 



= f o + —rr -I- Ps 

^ s=0 



(r — s)t 



(2.15) 



Here the first two terms are the free evolution of the test 
mass, the third (with the sum) is the influence of the 

meters' back-action forces (analog of radiation-pressure 
force in an interferometer), and the fourth, 

^r = - r f F{t')dt'dt = - r [rr - t')F{t')dt' , 
1^ Jo Jo A* Jo 

(2.16) 

is the effect of the classical force. The force F{t) 
is encoded in the sequence of classical displacements 
{(.i,£.2; ■■■■,^n}- It is also easy to show from the Heisen- 
berg equations that the meter's generalized coordinate 
after interaction with the test mass is 



Q''r'''°"^ -Xo-—rT-Y,Ps 
^ s=0 



(r — s)i 



(2.17) 



2. Vanishing of the output's self commutator 

The set of final meter coordinates Q = {Q^.Q^, 
Qjv-i} forms the final data string for data analysis. It 
has vanishing self commutator. 



[Qs, Qr] = for all s and 



(2.18) 



— a result that can be deduced from the vanishing single- 
measurement commutator [Xmeas ) Pafter] = [Eq. (2.8)] 
for the earlier of the two measurements. 

It is instructive to see explicitly how this vanish- 
ing commutator arises, without explicit reference to our 



single-measurement analysis. The test-mass contribu- 
tions to the Q's [xo and Po in Eq. (2.17)] produce 



[Qs^ Qr, 



test— mass — 



Po 

■X„ ST., -Xo 



Po 



i?i{r — s)t 



(2.19) 



which is the analog of Eq. (1.4) for an interferometer test 
mass. This must be cancelled by a contribution from the 
meters. Indeed it is. If (for concreteness) r > s, then 

the cancelling contribution comes from a commutator of 
(i) the 0°^°^°^"^ piece of Qs (the noise superposed on the 
output signal s by meter s) and (ii) the Pg term in Qr (the 
noise in the later measurement produced by the back- 
action of the earlier measurement): 



[Qs; Qr 



gbefore^ -P, 

—ih{r — s)t 



(r — s)t 



(2.20) 



In this example, one can trace these cancellations to 

the bilinear form xPg and xPr of each piece of the inter- 
action Hamiltonian. However, this type of cancellation 
is far more general than just bilinear Hamiltonians: In 
every sequence of measurements on any kind of system, 
by the time a human looks at the output data stream, 
its entries have all been amplified to classical size, and 
therefore they must all be classical quantities and must 
commute, [QsiQr] = [QsiQr] = 0. Remarkably, quan- 
tum mechanics is so constructed that, for a wide variety 
of measurements, the measured values (regarded as Her- 
mitian observables) commute even before the amplifica- 
tion to classical size. This is true in the above example. 
It is true in a realistic variant of this example involving 
pulsed- light measurements (Sec. II C). It is true in a vari- 
ant of this example involving continuous measurements 
by an electromagnetic wave in an idealized transmission 
line [32]. And, as we shall see in Sec. Ill A and Appendix 
C, it is also true for gravitational- wave interferometers — 
and indeed for all measurements in which the measured 
results are encoded in the photon number flux of a (quan- 
tized) electromagnetic wave; i.e., all measurements based 
on photodetection. More generally, it is true for any lin- 
ear measurement (Appendix B below, Ref. [3], and Eq. 
(2.34) of Ref. [17]); and, in fact, all the measurements 
discussed above, including gravitational-wave measure- 
ments, are linear. 

The classical nature of the output signal (the com- 
mutation of the data entries) guarantees that, when a 
human looks at one data entry, the resulting reduction 
of the state of the measured system cannot have any in- 
fluence on the observed values of the other data entries. 
Correspondingly, we can carry out any data processing 
procedures we wish on the Qr, without fear of introduc- 
ing new quantum noise. 
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3. Removal of test-mass influence from the output 

Our goal is to measure the classical force F{t) that 
acted on the test mass, without any contamination from 
the test mass's quantum properties — more specifically, 
without any contamination from uncertainty-principle 
aspects of the test mass's initial state. The initial state 
does influence the measured values Qr of the outgut ob- 
servables Qr, since in the Heisenberg Picture the Qr con- 
tain the test mass's initial position Xq and momentum Po 
[Eq. (2.f7)]. Therefore, our goal translates into finding 
a data analysis procedure that will remove from the out- 
put data set {Qi,Q2, ■ ■ •} all influence of the test-mass 
initial state (or equivalently all influence of Xq and Po), 
while retaining the influence of F{t). In fact, we can do 
so rather easily, regardless of what the test-mass initial 
state might have been. As we shall see, our ability to do 
so relies crucially on the linearity of our measurements; 
in particular, on the fact that the output observables Qr 
are linear in Xo and Po- 

To bring out the essence, we shall restrict ourselves to 
just three meters, N = 3. The generalization to large N 
is straightforward. 

The measured data sample Qr is equal to the freely 
evolving test-mass position at time rr, Xfree{t = rr) = 
Xo -\- (j)o/fJ-)rT (which is linear in Xg, Po), plus noise. 
Since the free evolution satisfies the equation of mo- 
tion cPxhfifi/dt^ = 0, it is a reasonable guess that we 
can remove the influence of Xo and Po from the data 
Qr by applying to them the discrete version of a sec- 
ond time derivative'^ (which is a linear signal process- 
ing procedure). Accordingly, from the measured values 
{Qo, Qi, Q2} of {Qo, Qi, Q2} in a representative experi- 
ment, we construct the discrete second time derivative 

R = {Q2- Qi) - (Qi - Qo) = Qo- 2Qi + Q2 • (2.21) 

The following argument shows that all the statistical 
properties of this quantity, in a large series of experi- 
ments (in which the initial states |in) of the test mass 
and meters are always the same) are, indeed, devoid of 
any influence of and Po, and thus are unaff'ected by 
the test-mass initial state.'* 



^In Sec. Ill C of Ref. [26], Caves uses his path-integral for- 
mulation of measurement theory to analyze measurements of 
the discrete second time derivative of the position of a free 
particle on which a classical force acts. His analysis reveals 
the same conclusion as we obtain in our pedagogical example: 
the measured quantity contains information about the force 
and is devoid of any influence from the particle's initial state. 

"^The crucial idea of avoiding the influence of the test- 
mass initial state by monitoring differences of observables 
[{Q2 — Qi) — {Qi — Qo) in our case] is contained in a pa- 
per and book by Alter and Yamamoto [33,34]. Alter and Ya- 
mamoto point out that, for a test mass on which a classical 



These statistical properties arc embodied in the means, 
over all the experiments, of arbitrary functions G{R). 
The theory of measurement tells us that, because the 
Q's all commute, the computed mean of G{R) is given 
by 



[computed mean of G{R)] = {m\G{R)\m) . 

where R is the operator corresponding to R 

R = Qo-2Qi + Q2 =-(^0-2^+6) + 

At 

I ^before Q^^hctore 1 _|_ ^before 



(2.22) 



(2.23) 



cf. Eq. (2.17). Because R is independent of Xo and Po, 

the computed mean (2.22) and thence all the measure- 
ment statistics of R will be completely independent of the 
test-mass quantum mechanics, and in particular indepen- 
dent of the test mass's initial state. Moreover, Eq. (2.22) 
implies that, so far as measurement results and statistics 
are concerned, measuring the Q's and then computing R 
is completely equivalent to measuring R directly. 
Although R is independent of Xo and Po it contains 



Co -2a +6 



M Jo 



T-\t-T\)F{t)dt = —F , 



(2.24) 



where is a weighted mean of the classical force F over 



force acts, the momentum p{t) at time t and the momentum 
p(0) at time are correlated in that p{t) = p{0) +/(,* dt'F{t'); 

so, if one measures p{t) — p{0) = dt' F{t'), one thereby can 
get information about the force without any contaminating 
influence of the test-mass initial state. They say (page 96 
of [34]) that this is so not only when one measures directly 

the difference pit) — p{Q) (as in Sec. 7.2.2 of their [34]), but 
also when the difference is determined computationally from 
the results of measurements of pit) and p(0) [an analog of 

our way of monitoring {Q2 — Qi) — (Qi — Qo)]- When going 
on to discuss position measurements, Alter and Yamamoto 

note that x{t) - x{0) = p{0)t/m -\- J* dt' J*' dt"F{t") /m, so 
a measurement of x{t) — x{Q) is contaminated [via p(0)t/m] 
by noise from the test-mass initial state. Examining this con- 
tamination, they conclude that "force detection via position 
monitoring of a free mass is limited by ... the SQL" [33]. 
While this conclusion is correct when one monitors x{t) — x{0) 
in the manner envisioned by Alter and Yamamoto, it is in- 
correct for the alternative strategy embodied in our model 
problem. Instead of monitoring x{t) — a;(0), one should moni- 
tor x{0) — 2x{t) -\- x{2t), which for a free mass is independent 
of both Xo = x{0) and po = p(0). Then the measurement 
output contains information about about the force F(t), un- 
contaminated by any influence of the test-mass initial state. 
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the time interval < t < 2r; cf. Eq. (2.16).^ Thus, this 
measurement of R is actually a measurement of F, and is 
contaminated by quantum noise from the meters but not 
by quantum noise from the test mass. The only role of 
the quantum mechanical test mass is to feed the classical 
signal F and the meter back-action noise Pir/m into the 
output. 

For those readers who are uncomfortable with our use 
of the Heisenberg picture to derive this very important 
result, we present a Schroedinger-picture derivation in 
Appendix A. 

This three-meter thought experiment is a prototype 
for our discussion of gravitational-wave interferometers 
in Sec. IIIB. There as here, the linearity of the output in 
the test-mass initial positions and momenta will enable us 
to find a linear signal processing procedure that removes 
the initial-state influence. Here that procedure was a dis- 
crete second time derivative. For an interferometer it will 
be a discrete Fourier transform of the measured photon 
flux (the output), and a discarding of Fourier compo- 
nents at the test masses' natural frequencies (the 1 Hz 
pendular swinging frequency in the case of conventional 
interferometers) . 

For an elegant path-integral analysis of the removal 
of test-mass initial conditions from the output of mea- 
surements of any harmonic oscillator on which a classical 
force acts, see the last portion of Sec. HI C of Caves [26]. 



4- The SQL for the classical-force measurement 

How small can the test-mass noise be? A "naive" op- 
timization of the meters leads to the standard quantum 
limit on the measured force, in the same way as a "naive" 
optimization of a gravitational- wave interferometer's de- 
sign (forcing it to retain the conventional LIGO-I opti- 
cal topology but optimizing its laser power) leads to the 
gravitational-wave SQL. Specifically: 

Let the three meters all be prepared in initial states 
that are "naive" in the sense that they have no correla- 
tions between their coordinates and momenta. Then Eqs. 
(2.23) and (2.24) imply that the variance of the measured 
mean force is 



(Ai^)^ = ^ (AQo^ 



-lbcforc\ 2 



(2AQ^'='"'=)2 



^Notice that, aside from meter noise, is equal to x{tr) — 
p(0)tr/n [Eq. (2.15)], which is a QND observable (as M.B. 
Mensky pointed out long ago). Therefore, the quantity II 
that we measure can be regarded as a discrete second time 
derivative of a QND observable which suggests that it can 
be the foundation for a QND measurement; see Sec. IIB5 
below. 



(2.25) 

Obviously, this variance is minimized by putting meters 
and 2 into (near) eigenstates of their coordinates, so 
^gbcforc ^ /^qY°''' = 0. To minimize the noise from 
meter 1, we require that it have the smallest variances 
compatible with its uncertainty relation, 



^gbefore^p^ = | , 



(2.26) 



and we adjust the ratio AQ\°^°'"^ / /S.Pi so as to minimize 
(AF)2. The result is 



{AFf = ^ 



(2.27) 



which is the SQL for measuring a classical force, up to a 
factor of order unity; cf. Sec. 8.1 of Ref. [3]. 

It is evident from this analysis that the true physi- 
cal origin of the SQL in classical force measurements is 
the meter's noise., not the test-mass noise. On the other 
hand, the quantum properties of the meter and of the 
test mass arc intimately coupled through the reqTiirement 
that the meter commutators cancel the test-mass commu- 
tator in the measurement output, so that [Qr,Qs] = 
[Eq. (2.18)]. This intimate coupling — which, as we have 
discussed, has enormous generality — ensures that the 
SQL can be derived equally well from test-mass consid- 
erations and from meter considerations. We saw this ex- 
plicitly in Sec. I for an interferometric gravitational-wave 
detector. 



5. Beating the SQL 

Equation (2.23) suggests a way to beat the classical- 
force SQL and, in fact, achieve arbitrarily high accuracy: 
As in our "naive" optimization, before the measurement 
we place meters and 2 in (near) eigenstates of their 
coordinates, so AQo = AQ2 = 0, but instead of putting 
meter 1 in a "naive" state with uncorrelated coordinate 
and momentum, we place it in a (near) eigenstate of 



^squeeze 



^before 



Ar/2/z . 



(2.28) 



(This meter- 1 state is analogous to the squeezed- vacuum 
state, which Unruh [18] has proposed be inserted into a 
conventional interferometer's dark port in order to beat 
the gravitational-wave SQL; see Sec. 11 C below.) These 
initial meter states, together with Eqs. (2.23) and (2.22). 
guarantee that the variance of the computed quantity R 
vanishes Ai? = 0, and thence [via Eqs. (2.24) and (2.23)] 
that the variance of the measured mean force vanishes, 
AF = 0. Thus, by putting the initial state of meter 1 into 
the analog of a squeezed vacuum state, we can achieve 
an arbitrarily accurate measurement of F. 

The SQL can also be evaded by modifying the meters' 
measured quantitites instead of modifying their initial 
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states. Specifically, measure Qq and Q2 as before, but on 
meter 1 instead of measuring the coordinate Qi, measure 
the following linear combination of the coordinate and 
momentum (with the coefficient a to be chosen below): 



gbefore _ _ PO ^ , 



T-aPi-^1. (2.29) 



From Eqs. (2.29), (2.17) and (2.18), we see that the out- 
put observables {Qo,QY''^ ,Q2} all commute with each 
other. Therefore, when we combine their measured val- 
ues into the discrete second time derivative 



2Qr + Q2 



(2.30) 



its statistics will be the same as if we had directly mea- 
sured the corresponding operator 



-Rvar = Qo — 2(3™'^ -|- Q2 

Qbefore _ 2gb 



-(Co -2a +6) 



+ 



before 2aPi + QY°''' 



(2.31) 



Evidently, we should choose 2a = r/fi, so the quantity 
measured is 



gr = Or + Pr^. 

Then Eqs. (2.31) and (2.24) imply that 



-Rvar — F + Qo 



before 



(2.33) 



Therefore, by measuring our chosen linear combination 
of meter 1 's coordinate and momentum, and then com- 
puting the discrete second time derivative, we have suc- 
ceeded in removing from our output observable i?var not 
only the test-mass variables x,,- po, but also the back- 
action influence of the meters on the measurement ( all 
three Pr 's)! Correspondingly, by putting the meters into 
"naive" initial states (states with no position-momentum 
correlations) that are near eigenstates of their coordi- 
nates (so AQo, AQi, A(52 are arbitrarily small and the 
back-action fluctuations APqi APi, AP2 are arbitrarily 
large), then from the computed quantity iZvar, we can in- 
fer the mean position F with arbitrarily good precision. 

This strategy was devised, in the context of optical 
measurements of test masses, by Vyatchanin, Matsko and 
Zubova [6-9], and is called a Quantum Variational Mea- 
surement. A gravitational- wave interferometer that uti- 
lizes it (and can beat the SQL) is called a Variational 
Output Interferometer [13]. 

Of course, one can also beat the SQL for force mea- 
surements by a combination of putting the meters into 
initially squeezed states and performing a quantum vari- 
ational measurement on their outputs. A gravitational- 
wave detector based on this mixed strategy is called 
a Squeezed VariaMonal Interferometer, and may have 
practical advantages over squeezed-input and variational- 
output interferometers [13]. 



C. Pulsed-light measurements of test-mciss position 

Our two pedagogical examples (single position mea- 
surement, Sec. II A, and classical force measm'cmcnt. Sec. 
II B) can be realized using pulsed-light measurements of 
the test-mass position. We exhibit this realization in 
part to lend reality to our highly idealized examples, 
and in part as a bridge from those simple examples to 
gravitational-wave interferometers with their far greater 
complexity (Sec. Ill below). 

In each pulsed-light measurement we reflect a laser 
light pulse, with carrier freriucncy lUq and Gaussian- 
profile duration Tq, off a mirror on the front face of the 
test mass, and from the light's phase change we deduce 
the test-mass position x averaged over the pulse. This 
is a concrete realization not only of the pulsed measure- 
ments of our pedagogical examples, but also of a Heisen- 
berg microscope. We presume that the pulse duration Tq 
is long compared to the light's period 27r/a>o, but short 
compared to the time r between measurements. 

We shall analyze in detail one such pulsed measure- 
ment. The electric field of the reflected wave, at some 
fiducial location, is 



(2.32) E{t) = 



2tt hojQ 



cS 



2iwn 



x{t) + a{t) 



h.c. , (2.34) 



where Aq is the pulse's amphtude, S is its cross sectional 
area, c is the speed of light, 2{LVo/c)x{t) is the phase shift 
induced by the test-mass displacement a;(t), "h.c." means 
Hermitian conjugate, and d{t) is the electric field's am- 
plitude operator. Because we are concerned only about 
timescales of order the pulse duration tq or longer, which 
means side-band frequencies ^ I/tq loq, wc can use the 
quasimonochromatic approximation to the commutation 
relation for a{t) [35]: 



[a{t),a\t')] =5{t-t') . 



(2.35) 



Note that, when decomposed into quadratures with re- 
spect to the carrier frequency, this electric field is 



E{t) = EA{t) cos LJot + E^{t) sin Wot , 



(2.36) 



where Ea and E^, the amplitude and phase quadratures 
(i.e., the quadrature components oriented along and per- 
pendicular to the amplitude direction in the quadrature 
plane) are given by 



Ea 

Ea, 



= 2 



2-KhUo 

cS 

2'Kh,UJo 

cS 



Aoe~ 



eiirl , ( a(t) +aU^) 



2Ao—e 

c 



+ 



a{t) - a) it) 
2i 



, (2.37a) 



(2.37b) 
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The power W{t) in the incident wave can be written 
as the sum of a mean power {W{t)) and a fluctuating 
(noise) part W{t): 



[^meas ? Pafter] — 



(2.44) 



Wit) = Sc^^ = (Wit)) + W{t), 



47r 



{W{t))=hojoAle-''/^o^ 
W{t) = 2nujQ Ao e-*'/=^"'o 



(2.38a) 
(2.38b) 
(2.38c) 



Here the over bar means "average over the carrier pe- 
riod". The hght-pressure force on the mirror is F{t) = 
2W{t)/c. The fluctuating part of this, F{t) = 2W(t)/c, 
is the back-action of the measurement on the test mass, 
and it produces the back-action momentum change 



BA 



/- 



dt 



2W{t) 



oo 

47ia;i 



The test-mass momentum before and after the pulsed 
measurement are related by 



Pafter = Pbefore + ^PBA • 



(2.40) 



The experimenter deduces the phase shift {2LUo/c)x{t) 
and thence the test-mass displacement x{t) by measuring 
the electric field's phase quadrature (e.g., via inter- 
fcrometry or homodync detection). More precisely, the 
experimenter measures the phase quadrature integrated 
over the pulse, obtaining a result proportional to 



cS 



2hWo 4:'JTUoToA{ 



/-\-oc 
-oo 

(2.41) 



cf. Eq. (2.37b). Here x is the mirror position averaged 
over the short pulse, fmeas is the measured value of x, 
and Sxmeas is the measurement noise superposed on the 
output by the light pulse 



dXr, 



2^/^^uJoTQAo \ 2i 



(2.42) 



It is straightforward, from the commutator \a{t), a^{t')] = 
d{t ~ t'), to show that the measurement noise and the 
back-action impulse have the same commutator 



[^^measj^PBA] = -ih 



(2.43) 



as for the idealized single measurement of Sec. II A [Eq. 
(2.7)], and correspondingly the mirror's measured posi- 
tion and its final momentum commute. 



The fundamental equations (2.41), (2.40), (2.43) and 
(2.44) for this pulsed-light measurement are the same as 
those (2.4), (2.7), (2.8) for our idealized single measure- 
ment, and this measurement is thus a realistic variant of 
the idealized one. Similarly, a sequence of pulsed-light 
measurements can be used to monitor a classical force 
acting on a mirror, and the fundamental equations for 
such measurements are the same as for the idealized ex- 
ample of Sec. II B. 

In such pulsed-light experiments, the measurement 
noise dxmc&s is proportional to the fluctuations of the 
light's phase quadrature [Eqs. (2.37b) and (2.42)], 
and the back-action impulse Spba is proportional to the 
fluctuations of its amplitude quadrature Ea [Eqs. (2.37a) 
and (2.39)]. Of course, experimenters can measure any 
quadrature of the reflected light pulse that they wish. 
To achieve a QND quantum variational measurement of 
a classical force acting on the test mass [6-9], the ex- 
perimenter should measure Q™'' = Qi + P\t/2ijl in the 
language of our idealized thought experiment [Eq. (2.32)], 
which [by Eqs. (2.5)] translates into — (S.inicas + <^PBAT/2/i 
plus the light's signal and carrier, which in turn is a spe- 
cific linear combination of the light's amplitude and phase 
quadratures Ea and E^ [Eqs. (2.37), (2.42), (2.39)]. The 
experimenter can also prepare the incident pulse in a 
squeezed state, in the manner required for an Unruh- 
type [18] QND measurement of the classical force. In 
the language of our idealized thought experiment, the de- 
sired squeezed state is a (near) eigenstate of Q^i"*^*^^** = 

Qi — Pit/2ijl [Eq. (2.28)], which translates into a near 
eigenstate of (5xmeas + <5pBAT/2/x [cf. Eqs. (2.5)], or equiv- 
alently a near eigenstate of a specific linear combination 
of Ea and Ea,. 



III. GRAVITATIONAL-WAVE 
INTERFEROMETERS AND OTHER 
PHOTODETECTION-BASED DEVICES 

We now turn our attention to gravitational-wave in- 
terferometers and other real, high-precision devices for 
monitoring classical forces that act on test masses. Our 
goal is to prove that for these devices, as for our lAc- 
alized examples, the force-measurement precision can be 
made completely independent of the test mass's quantum 
properties, including its initial state and that this can be 
achieved by an appropriate filtering of the output data 
stream. 

As in our examples, this conclusion relies on the van- 
ishing commutator of the observables that constitute the 
output data stream. We shall now discuss the nature of 
the output data stream and show that its commutator 
does, indeed, vanish. 
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A. Vanishing commutator of the output 

For interferometers and many other force-monitoring 

devices, the data stream, shortly before amphfication to 
classical size, is encoded in an output light beam, and 
that beam is sent into a photodetector which monitors 
its photon number flux Af(t) . The photodetector and as- 
sociated electronics integrate up /^{t) over time intervals 
with duration r long compared to the light beam's car- 
rier period, r ^ 2tt/uJo ^ 10^^^ s, but short compared 
to the shortest timcscalcs on which the classical force 
changes (r <C tqw ~ 10~^ s for the gravitational waves 
sought by interferometers). For LIGO-I interferometers, 
the integration time has been chosen to be r = 5 x 10~^ 
s. The result is a discrctized output data stream, whose 
Hermitian observables are the numbers of photons in the 
successive data samples. 



-f 

O — ( 



s{t-tj)Af{t)dt . 



(3.1) 



Here tj ~ jr^ is the time of sample j, and s{t) is a sam- 
pling function approximately equal to unity during a time 
interval At = Tq centered on tj and zero outside that time 
interval. 

The photon number samples Nj are the analogs, for an 
interferometer or other force-monitoring device, of the 
meter coordinates Qj in the idealized example of Sec. 
II B. 

In Appendix C we show that for any free light beam, 
the number flux operator, evaluated at a fixed plane or- 
thogonal to the optic axis (e.g. at the entrance to the pho- 
todetector) self commutes, 



m),j^{t')] = o. 



(3.2) 



This guarantees, in turn, that all the output photon- 
number data samples (3.1) commute with each other 



[Nj,Nk]=0. 



(3.3) 



As we shall see below [Eq. (3.9)], the initial position and 
momentum of the test mass, Xo and po, appear linearly 
in the output variables Sf{t) and Nj. They obviously will 
produce nonzero contributions to the output commuta- 
tors. As in our simple examples (Sec. II), these nonzero 
test-mass contributions must be cancelled by identical 
nonzero contributions from noncommutation of the mea- 
surement noise (photon shot noise) and the back-action 
noise (radiation- pressure noise). 



B. Devising a filter to remove test-mass quantum 
noise 

The vanishing output commutators constitute our first 
underpinning for freeing the measurements from the 



influence of test-mass quantization. As in the ideal- 
ized measurements of Sec. II B, the vanishing commu- 
tators guarantee a key property of the data analysis: If, 
from each specific realization of the output data stream 
{A^i, A^2, • . •}, our data analysis produces a new set of 
quantities (the "filtered output variables") 



RjiNr,N2,...) 



(3.4) 



then the statistics of these Rj will be identically the same 
as if we had directly measured the corresponding observ- 
ables 



Rj{Ni,N2,...) 



(3.5) 



rather than computing them from the measured Nj's. 

Therefore, we can regard our interferometer (or other 
device) as measuring the filtered output observables 
i?2, . . .}, whatever those observables may be. 
By analyzing the test-mass dynamics of the inter- 
ferometer (or other measuring device) in the Heisen- 
berg picture, one can learn how the test-mass initial 
position Xo and momentum po influence the operators 
{Ni, N2, ■ ■ ■}■ One can then deduce a set of filtered ob- 
servables {Ri,R2, . . .} in which Xo and Po do not ap- 
pear but the gravitational-wave or other classical force 
information is retained. (These will be the analogs of 
R = - 2Qi + Q2 [Eq. (2.23)] in our simple model 
problem). The filter that leads from {A^i, A^2, • • •} to 
{Ri, R2, . . .}, when applied to the output (c-number) data 
{Ni,N2, ■ ■ ■} to produce {Ri,R2, ■ . .}, is guaranteed to 
remove all influence of Xq and po, and thence all influ- 
ence of the test-mass initial state. 



1. Influence of Xo and po on the output data 

To make this more specific, let us explore how Xq and 
Po influence the output data train. 

To very high accuracy (sufficient for our purposes), in- 
terferometers (and most other force-measuring devices) 
are linear. The inputs are: (i) the test-mass position x{t) 
[actually, the difference between four test-mass positions 
in the case of an interferometer; Eq. (1.2)], and (ii) the 
electric field operators Ea{t), a = 1,2,... for the field 
fluctuations that enter the interferometer at the bright 
port, at the dark port, and at all light-dissipation loca- 
tions (e.g., at mirrors where bits of light scatter out of 
the optical train and reciprocally new bits of field fluctu- 
ations scatter into it); see, e.g., the detailed analysis of 
interferometers in Ref. [13]. The output photon flux is a 
linear functional of these inputs, 



m = f 

J —I 



K,{t - t')x{t') + - t')Ea(t') 



dt' ; 
(3.6) 
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cf. the discussion in Appendix C. The Ea terms con- 
stitute the photon shot noise (analogs of Q^^^™"^ in our 
ideahzed example, Sec. II B). 

The test-mass initial obscrvablcs Xg and po enter Af{t) 
and thence {ATi, iV2, . . .}, through x(t) in a manner gov- 
erned by the test masses' free dynamics. The nature of 
that free dynamics depends on the interferometer design. 
We shall consider two examples in turn: interferometers 
with pendular dynamics, and signal-recycled interferom- 
eters. These examples should be easily extendable to any 
other type of interferometer than might be conceived in 
the future. 



2. Interferometers with pendular dynamics 

In conventional gravitational-wave interferometers 
(e.g. LIGO-I, VIRGO and TAMA) and in the QND inter- 
ferometers analyzed by Kimble et. al. [13], the test masses 
swing sinusoidally at ~ 1 Hz frequency in response to 
their suspensions' pendular restoring force (as modified 
slightly by the optical cavities' radiation-pressure force): 



Xfieeit) = Xo COS Umt " 



Po 



(3.7) 



Here fi is the reduced mass (1/4 the actual mass of one 
test mass in the case of an interferometer) and utm ~ 
27r X 1 Hz is the pendular swinging frequency. There is no 
significant damping of the free motion (3.7) because the 
experimenters take great pains to liberate the test masses 
from all damping; the typical damping times in LIGO- 
I are of order a day, and in advanced interferometers 
(LIGO-II and beyond) will be of order a year or more 
[14,36], which is far longer than the data segments used 
in the data analysis. 

Superimposed on the free test-mass dynamics (3.7) arc 

(i) the influence ^gw {t) of the gravitational- wave signal, 

(ii) the "back-action" influence XBA{t) of the light's fluc- 
tuating radiation pressure (which is linear in the input 
fields Ea and is the analog of the Pr and SpsA of our dis- 
crete model problems), and (iii) the infiuence Mother (0 of 
a variety of other forces — low-frequency feedback forces 
from servo systems, thermal-noise forces, seismic vibra- 
tion forces, etc: 

X{t) = ffree(t) + ^Gw{t) + XBA{t) + UheAt) • (3.8) 

Inserting Eq. (3.7) into (3.8) and then (3.8) into (3.6) 
we see that, for a test-mass with pendular dynamics, the 
initial test-mass position and momentum operators ap- 
pear in the output flux operator in the form 



Af{t)= f K^t-t') 



Xo cos Ulmt' H — sin LJmt' 



-|- (other contributions) . 



dt' 
(3.9) 



The interferometer's transfer function Kx{t — t') is inde- 
pendent of absolute time and thus transforms frequency- 
ojm inputs into frequency-o.)™ outputs. Therefore, Xo and 
Po appear in the output solely at frequency u>ml'2''!^ ~ 1 
Hz. Now, because the output data generally have large 
noise (seismic and other) at frequencies below ~ 10 Hz, 
it is routine, in interferometers, to high-pass filter the 
output data so as to remove frequencies below ~ 10 Hz. 
When one does so, one automatically removes all influ- 
ence of Xo and Po from the filtered data Rj [Eq. (3.4)]. 
This is a precise analog of applying the discrete second 
time derivative to the output data in our simple exam- 
ple (Sec. II B) so as to remove Xo and Po from the data; 
and it is a realization of a general class of measurement 
procedures, for a harmonic oscillator on which a classi- 
cal force acts, that is analyzed by Caves using his path 
integral formalism (last part of Sec. Ill C of Ref. [26]). 



3. Signal-recycled interferometers 

A signal-recycling mirror, placed at an interferometer's 
output port, sends information about the test-mass po- 
sition x{t) back into the interferometer as part of the 
back-action (radiation-pressure) force, and thereby alters 
the free test-mass dynamics. The altered free dynamics 
have been analyzed in detail by Buonanno and Chen [17]; 
they find that the test masses and the interferometer's 
side-band light form a coupled system with four degrees 
of freedom, so Xq and po appear in XfT-cc{t), and thence 
in x{t) and thence in J\f{t) at four discrete frequencies 
(jjA {A = 1,2,3,4). Correspondingly, in the output data 
train, the influence of the test-mass initial state is con- 
fined to the Fourier components at the frequencies lja- 

If these frequencies were real, then one could remove 
the influence of the test-mass initial state from the data 
by filtering out the data's Fourier components at these 
four frequencies. However, as Buonanno and Chen [17] 
discuss, such filtering is not necessary: The frequencies 
arc actually complex with imaginary parts that produce 
damping on timescales ;$ 1 second (when a servo is intro- 
duced to control an instability). Therefore, the influence 
of Xo and po on the output flux operator Sf{t) damps out 
quickly, and correspondingly (see the end of Sec. Ill A), 
the influence of the test-mass initial state on the output 
data train damps out quickly without any filtering. 



IV. CONCLUSIONS 

To reiterate: In an interferometer (and many other 
force- measuring devices), the output signal is encoded in 
the photon number flux operator M{t) of a light beam, 
which is converted into discrete photon number sam- 
ples Nj by a photodetector and electronics. These out- 
puts have vanishing commutators [N{t),N{t')] = and 
[Nj,Nk\ = and thus can be thought of as classical 
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quantities. These outputs arc linear in the initial test- 
mass position Xo and momentum po and involve no other 
test-mass variables. The output commutators manage to 
vanish because the photon back-action noise and photon 
shot noise have commutators that cancel those of Xo and 

Po- 

In the output Af{t) of any interferometer with pen- 
dular dynamics, Xo and po appear only at the pendular 
frequency ujm/2n ^ 1 Hz, and all influences of Xq and po 
(including all influences of the test-mass initial state) are 
removed completely from the data by the high-pass filter- 
ing that is routine for interferometers. For other types of 
interferometers, with different test-mass dynamics, other 
data filtering procedures will remove the influence of Xg 
and Po and the test-mass initial state — and in some 
cases (e.g., a signal-recycled interferometer) no filtering 
is needed at all. 

This complete removal of all influence of Xo and po from 
the filtered data implies the answers to the three ques- 
tions posed in the introduction of this paper (Sec. I): 
(i) The test-mass quantum mechanics has no influence 
on the interferometer's noise; the only quantum noise is 
that arising from the light, (ii) Therefore, when ana- 
lyzing a candidate interferometer design, one need not 
worry about the test-mass quantum mechanics, except 
for using it to feed the gravity-wave signal and the back- 
action noise through the test mass to the photon-flux 
output, (iii) Similarly, when conceiving new designs for 
interferometers, one need not worry about the test-mass 
quantum mechanics — except for devising appropriate 
data filters to remove Xo and Po from the data. 
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APPENDIX A: TRIPLE MEASUREMENT IN 
THE SCHROEDINGER PICTURE 



force F acting on the test mass without any contaminat- 
ing noise whatsoever from the test mass's initial state. 
Our analysis will proceed in three steps: (i) an analy- 
sis of one of the position measurements (any one of the 
three). Sec. Al; (ii) [relying on step (i)] a derivation of 
the probability density W{Qo,Qi,Q2) for the outcome 
of the triple measurement procedure, Sec. A 2; and (iii) a 
use of this probability density to show that the combina- 
tion R = Qo — 2Qi + Q2 of the measurement results con- 
tains the desired information about F uncontaminated 
by any noise from the test-mass initial state, Sec. A3. 



1. Single position measurement 

Let |\E') be the state of the test mass before the mea- 
surement and 



f 

J — C 



mm dQ 



(Al) 



be the initial state of the meter, where the meter's eigen- 
states are normalized by 



{Q'\Q) = S{Q - Q') 



(A2) 



We leave the test-mass state l^*) completely unspecified 
since our goal is to show that it has no infiuence at all on 
the measurement outcome. For concreteness we specify 
the meter's initial wave function V'(<3) to be Gaussian: 



exp 



n ) 



(A3) 



Here Ag (denoted AQ''<'f°''<' in the text) is the initial 
variance of Q and 



Aqp = 



(gP + PQ) 



(A4) 



is the initial cross correlation of the meter's position and 
momentum. For this Gaussian initial state, the variance 
Ap of the meter's momentum (denoted ^p^'^^°^'^ in the 
text) is given by the minimum-uncertainty relation 



A 2 A 2 A 2 



(A5) 



The first stage of the measurement process is the 
interaction of the test mass and the meter. In the 
Schroedinger Picture this interaction puts the meter and 
test mass into the entangled state 



In this appendix we present a Schroedinger-picturc 
analysis of the most important of this paper's pedagogical 
thought experiments (Sec. II B): a triple measurement of 
the position of a free test mass, using three independent 
meters, with the goal of determining the mean classical 



where 




(A6) 



(A7) 
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is the evolution operator associated with the interaction 
(delta function) part of the Hamiltonian (2.14). 

The next stage is a precise measurement of the meter's 
generalized position Q. This measurement disentangles 
the quantum states of the test mass and meter: the meter 
gets reduced to the eigenstatc \Q) of Q, where Q is the c- 
number obtained as result of this measurement, and the 
test mass gets reduced to the state 



{Q\um^) mm 



W{Q) 



W{Q) 



where 



n(Q) = (Q|[7|V) 



(A8) 



(A9) 



is the reduction operator describing the entire two-stage 
measurement procedure, and 



W{Q) = mh\Q)h{Q)\^) 



(AlO) 



is the probability density for obtaining the result Q. 

An explicit form for the reduction operator can be ob- 
tained by substituting Eqs. (Al), (A3) and (A7) into Eq. 
(A9); the result is: 



^{Q) = (Qlexp 



ixP 



mm dQ 



= {Q\ / \x){x\m)\Q-x)dxdQ 

J — oo 

/oo 
\x){x\'il){Q + x) dx 
-oo 



1 



exp 



Q 



{Q + xTfl 
2A| ^2 



n 



(All) 



where we have used the shift-operator relation 
e^^/^'IQ) = \Q-x)= J^^ dx\x){x\Q - x) and the rela- 
tion (QIO -x) =5{Q-x~ Q). 

We will need below the following formulae (some are 
evident, and for the others we provide outlines of the 
proofs): 



/oo 
fit(Q)n(Q)dQ: 
-OO 



1 



h'<{Q)(l{Q)QdQ = -x 

) 

/OO 
-OO 



x^+A^ 



Q ' 



(A12) 
(A13) 

(A14) 



/OO 
nt(Q)i"n(Q)rfQ = i" 
-oo 



(n = 0, 1, . . .) , (A15) 



/OO 
il'' {Q)xn{Q) Q dQ 
-OO 



(A16) 



n\Q)pn{Q) dQ 

n\Q) (n{Q)p+ [p,f2((3)]) dQ 



= r n\Q)n{Q)dQp-ih r h{Q)^^^^j^dQ 



= p, 



(A17) 



/OO 
n^{Q)p''Cl{Q)dQ 
-oo 

= f to(Q) + \n\Q),p\) (Cl{Q)p+ \p,n{Q)]) dQ 

J —oo 

/OO pc 
h\Q)h{Q)dQp + 'h' / 
-oo J — c 



dCt^iQ) dn{Q) 



dQ 



1 fn^ 



/oo 
h\Q)ph{Q)QdQ 
-oo 

= f n\Q) (n{Q)p+ \pM{Q)\) QdQ 

J — oo ^ L J/ 

= r Cl\Q)CliQ)QdQp-th r Cl(Q)^^}l^QdQ 

J —oo J — oo 

... .^/l Aqp\ xp+px 
= -xp + ihi- — ) = 2 — + ^"3^' 



/oo 
n\Q){xp + px)Cl{Q) dQ = xp- 
-oo 



px 



(A20) 



2. The triple measurement procedure 

The triple measurement procedure described in Sec. 
II B 1 of the text consists of the following five stages. 

1. An initial position measurement of the type we have 

just analyzed, using meter number 0. This mea- 
surement reduces the test mass's wave function to 



f^o(Oo)|*) 



WoiQo) 



(A21) 



[Eq. (A8)], where Oo(Qo) is the reduction operator 
[Eq. (A9)], and Qo is the result of this measure- 
ment. The probability density for obtaining this 
result is equal to 
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WoiQo) = (*|OS(Qo)Oo(go)|*) (A22) 
[Eq. (AlO)]. 

2. Free evolution of the test mass during the time r. 

Denoting the corresponding evolution operator by 
Uq, the test-mass wave function after this stage is 
given by 



WoiQo) 



(A23) 



Equation (A28) is the principal result of this subsec- 
tion. We shall use it to study the statistics of the mea- 
surement outcomes. In that study we shall need the fol- 
lowing expression for each of the three reduction opera- 
tors [Eq. (All)]: 



1 



exp 



27rA 



Qs 



where s = 1, 2, 3. 



(Qs + 



QPs 



(A29) 



3. Second position measurement of the same type as 
in the first stage, but using a new meter, number 
1. The measurement result is denoted Qi, the re- 
duction operator is and the measurement 
reduces the test-mass state to 



Oi(Qi)40o(Qo)|^) 
^VFi(Qo,Qi) 



(A24) 



where 



WiiQo,Qi) 

= {^\Cll{Qo)U^onl{Qi)ni{Qi)UoMQo)\'^} (A25) 

is the joint probability disrtibution for the first two 
measurement results, Qo and Qi. 

4. Second free evolution of the test mass with the evo- 
lution operator Ui. After this stage the test-mass 
wave function is 



Z^iOi(QiHOo(Qo)|^) 
^Jw^{Qo,Qi) 



(A26) 



5. Finally, a third position measurement using a new 
meter, number 2, with the result Q2- After this 
measurement the test-mass state is 



fi2(02)Wifii(Qi)Wofio(Qo)|*) 



W2{Qo,Ql,Q2) 



(A27) 



where 



W2{Qo,Qi,Q2) = {■^\nl{QoM^\{Qi)ulnl{Q2) 

xd2{Q2)UMQi)UoMQo)\^) ■ (A28) 

is the joint probability distribution for all three 
measurement outcomes. 



3. Statistics of the measurement results 

If an explicit form for the initial wave function 
were specified, then the probability density (A28) could 
be calculated directly. However, that calculation would 
be very cumbersome, the final result would be quite com- 
plicated, and we have no need for it. Our final goal is 
not to study W2, but rather to analyze the statistics 
of the quantity R = Qo — 2Qi + Q2, which the exper- 
imenter computes from the three measurement outcomes 
Qs after the triple measurement procedure is complete. 
Specifically, we wish to verify the results of the text's 
Heisenberg-picture analysis: (i) That the mean value of R 
over a large number of experiments is (R) = (— T^//i)F, 
where r is the time between each pair of measurements, 
fi is the mass of the test mass, and F is the mean force 
that acts on the test mass [Eqs. (2.24) and (2.23) of the 
text], (ii) That the variance of R (and thence of the mea- 
sured value of F) is independent of the test-mass initial 
state I*), and is given by Eq. (2.25) when the meters' 
individual initial states have no position-momentum cor- 
relations, Aqps = 0, and can be made to vanish by a 
clever, "squeezed" choice of the meters' initial states. 

a. Mean value. The mean value of R over a large 
number of experiments is determined by the joint prob- 
ability distribution VF3 for the measurement outcomes: 

(i?) = (Qo - 2Qi + Q2) 

/OO 
(Qo - 2Qi + Q2)W2(Qo, Qi, Q2) c/QorfQidQa- 
-00 

(A30) 

Using Eqs. (A12), (A13), we bring this into the form 

/OO 
{^\nl{QoMn\iQiPl 
-00 

X (Qo - 2Qi - Wifii(Qi)Z^ofio(Qo)|*) dQodQi . 

(A31) 

Taking into account that 



UlUi 



1, 



pT 



U]xU\ = x+- \-xfi , 



(A32) 
(A33) 
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where /U is the mass of the test mass and 



1 f^^ 

XFi^ - / {2t -t)F{t)dt 
M Jt 



(A34) 



is the displacement of the test mass during stage 4 (the 
second interval of free evolution) caused by the external 
force F{t), expression (A31) can be further reduced to 
the form 

/oo 
-oo 

-x-^-XFi) Oi(Qi)Z^of^o(Qo)|*) dQodQi . 
M / 

(A35) 

The next calculations are just a repetition of the previous 
ones, with only the addition of Eqs. (A15), (A17) and 



pT 

U^pUo =P + Pfo 



t PT 

UqxUo = x-\ h a;_Fo , 



(A36) 
(A37) 



where 



xfo = - [ (r - t)F{t) dt , (A38) 
Jo 

Pfo= [ F{t)dt . (A39) 
Jo 

They give: 

{R) = J {^\nl{Qo)Ul (qo + 2x-x-^-xfi^ 



-L 



xUodoiQo)]^) dQo 

OO 

■oo 

X { Qo + X + xfo 



PFOJ 



-XFi]no{Qo)\^)dQo 



/Til PFOT \ I T \ Pfot 
(w Xfo XFi w) = XFO xfi 



■- 



{T-\t-T\)F{t)dt = F. (A40) 



This agrees with the Heisenberg-picture prediction [Eqs. 
(2.24) and (2.23) of the text, where we must note that 
the meters' initial states have {Qs) = (P?) = 0]. 

b. Variance The mean square value of the measure- 
ment outcome R over a large number of experiments is 
given by 

{R") = {{Qo - 2Qi + Q2f) 

/OO 
(Qo - 2Qi + Q2fW2{Qo, Qi, (32) dQodQidQ2 . 
-OO 

(A41) 



Using Eqs. (A12)-(A20), (A32), (A33), (A36), and 
(A37), we obtain: 

{R^) = ((Q0-2Q1 + Q2)') 

/OO 
{<i/\nl{QoMnl{Q,)Ul {Qo - 2Qi - xf + A| 
-00 

xUiQi{Qi)UoClo{Qo)\^) dQodQi 

/OO 
(*|f)S(Oo)Wo^f^l(Qi) 
-00 

pr 



Qo - 2Qi - x - Xfi 

m 



xOi(Qi)Zioi^o(Oo)|*) dQodQi 

/OO 
{mliQoP^o 
-00 



Qo + x Xfi 

m 



+ 4A2 



m 

00 



{ml{Qo) 

4AQpir 



Qo + -i + XF 
Apir^ ' 



Z^o^o(Qo)|*) ciQo 
Pfot 



-Xfi. 



fio(Qo)|*)rfQo 



{xfo 



PfoT 



m 

4Aqp it 



ApiT 



m \ M 

2 , a2 



+ A2 



Q2 



I*) 



= (Qo-2Qi + Q2)' + A^o 



+4A 



4AQPir 



m 



+ 



AP1^ \ , a2 

+Aq2 



(A42) 



Subtracting off the square of the mean, (i?)^ = 
(Qo — 2Qi + (52)^, we obtain for the variance of the com- 
puted quantity R, over many experiments. 



,(AF)2 = {ARf = {R^) - {Ry 



= A|o+4A^i 



4A, 



QPlT 



m 



+ 



ApiT 



a: 



2 1 

(A43) 



see Eq. (A40) for the first equality. This variance is inde- 
pendent of the test-mass initial state |^), in accord with 
prediction of the Heisenberg-picture analysis [passage fol- 
lowing Eq. (2.22) of the text]. When the three meters 
are all prepared in "naive" initial states, i.e. in states 
with imcorrelated generalized position Qs and momen- 
tum Pg, i.e. when Aqp,, = 0, then the variance (A43) has 
the form that we deduced using the Heisenberg picture 
[Eq. (2.25) ]. When the meters are prepared in the more 
clever "squeezed" manner, i.e. in near eigenstates of Qo, 
^squeeze = _ p^T-/2/i and Q2, then the variance (A43) 
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vanishes, in accord with the Heisenberg-picture predic- 
tion [passage following Eq. (2.28)]. 



APPENDIX B: LINEAR MEASUREMENTS 

An important feature of onr pedagogical examples 
(Sec. II), and of measurements performed by interfero- 
metric gravitational-wave detectors, is that they all are 
linear measurements in the sense of Ref. [3]; i.e., they all 
satisfy the following two conditions: 

(i) Linearity of the output: The meter's output can be 
written as the sum of the operator for the test object's 
measured variable and the operator for the meter's addi- 
tive noise [cf. Eq. (2.4a)], and the additive noise does not 
depend on the initial state of the test object. Formally 
this sum is an operator, but it can be treated as a classi- 
cal variable because it turns out to commute with itself 
at different times. 

(ii) Linearity of the back action: The measurement- 
induced perturbations of all the test-object observables 
that are involved in the measurement procedure can be 
described by linear formulas similar to Eq. (2.4b), and 
the perturbations [e.g. the second term on the right side 
of (2.4b)] do not depend on the initial state of the test 
object. 

This second condition requires discussion: The per- 
turbations' independence of the test-object initial state 
is particularly important when several test-object vari- 
ables are measured consecutively for example, if the 
same Heisenberg-Picture variable is measured quickly 
and repetitively at different moments of time as in our 
pedagogical examples (Sec. II), or if a variable is mea- 
sured continuously as in a gravitational-wave detector 
(Sec. III). Suppose, for example, that the variable xi is 
measured with precision Ax^^^ thereby perturbing, via 
back-action, some other variable £2- Then the accuracy 
of a subsequent measurement of £2 will be constrained 
by the perturbation 



Ax^' 



2Aa;f 



rl([^i'^2])|. 



(Bl) 



Our condition (ii) of back-action linearity requires that 

this perturbation not depend on the initial state of the 
test object. A sufficient condition for this is that the 
commutator [a;i,a;2] be a c-number, and that this re- 
quirement be fulfilled for all the operators involved in 
the measurement® 

Linear measurements are closely related to linear sys- 
tems (those for which the equations of motion for the 



®It can be shown that a slightly weaker condition is suf- 
ficient: second-order commutation of all these operators, 
[xi, [xj,Xk]] = for all i,j, k. 



generalized coordinates and momenta are linear; for ex- 
ample, a free mass and a harmonic oscillator) because the 
commutators of such systems' coordinates and momenta 
are c-numbers. 

In nonlinear measurements (e.g. measurements of a 
particle in a double- welled potential), some very strange 
phenomena can arise, for example the quantum Zeno ef- 
fect. 

Strictly speaking, all real meters are nonlinear. How- 
ever, in most cases they can be regarded as linear to high 
accuracy. For example, if one measures displacements 
of a mirror of a Fabry-Perot cavity by monitoring the 
phase of light that passes through the cavity (as is done 
in LIGO), then the measurements are linear so long as 
the displacements are much smaller than the width of a 
cavity resonance, i.e. much smaller than X/J^ where A is 
the wavelength of the light and !F is the cavity finesse. 

If, by contrast, the displacements are comparable to 
or much larger than X/!F, then the measurements are 
strongly nonlinear. An example is a proposed null- 
detector technique [37] for measuring the phase of a me- 
chanical oscillator, in which the oscillating mass is an end 
mirror of a Fabry-Perot cavity, and the times at which 
the mirror passes through cavity-resonant positions are 
measured with high accuracy by the cavity's momentary 
transmissivity. These measurements are highly nonlin- 
ear because, in the proposed design, not only are the 
mirror displacements large compared to the cavity's lin- 
earity regime, X/J^; the mechanical oscillator's amplitude 
of zero-point oscillations Sx^p is also large compared to 
A/JF. State reduction plays an important role in this 
null detector's measurements: it drives the mechanical 
oscillator into a squeezed-phase state, thereby facilitat- 
ing a high-precision monitoring of the oscillator's phase 
[37]. It would be instructive to analyze the use of this 
highly nonlinear meter to monitor a classical force that 
acts on the oscillator's mass. Does the oscillator's initial 
quantum state influence the accuracy of the monitoring? 

Three properties of an interferometric gravitational- 
wave detector (interferometric position meter) allow one 
to consider it as linear with sufficiently high precision to 
justify the linear analysis given in this paper. First, its 
test-mass mirrors can be regarded as free masses (or as 
harmonic oscillators if significant electromagnetic rigid- 
ity exists in the system). Second, its linearity range 
X/J^ ^ lO^^cm is much greater than the wave-induced 
displacements of the test masses (< 10~^^ cm). Hence, 
the signal phase shift of the output optical beam depends 
linearly on the displacement. Third, the measurement of 
the photon flux out the dark port is virtually equivalent 
to the measurement of the phase of the output beam 
because (i) the signal phase shift is much less than one 
radian and (ii) the mean value of the amplitude of the 
optical pumping field is much larger than the quantum 
uncertainties of its quadrature amplitudes. 

For a detailed presentation of the theory of linear 
measurements see Chaps. 5 and 6 of Ref. [3]. For 
a detailed application of this theory to interferometric 
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gravitational- wave detectors see Ref. [17]. 



does not self-commute, 



APPENDIX C: VANISHING SELF 
COMMUTATOR OF THE PHOTON NUMBER 
FLUX 

For any light beam (or other electromagnetic wave with 
confined cross section) , the number flux operator at some 
chosen transverse plane (e.g. the entry to a photodetec- 
tor) is 

A^W=/°°^/°"^4a.'e^^— (CI) 

Here al^ is the creation operator and au, the annihilation 
operator for photons of frequency oj, and their commu- 
tators are 

a^'] = [al, al,] = , [oo,, aj,,] = 27r5(w - w') . 

(C2) 

It is straightforward to veriiy from Eqs. (CI) and (C2) 
that 

[Mt),mt')]=0. (C3) 

Although this result is completely general, it is instruc- 
tive to derive the vanishing self commutator for the spe- 
cialized type of light beam that is used in interferometers 
and other force- measuring devices: a beam consisting of a 
monochromatic carrier with frequency plus sidebands 
embodied in a^. and aj^ . In this case to high accuracy we 
can linearize in the product of the carrier field and the 
side-band fields, obtaining for the relevant (side-band) 
photon flux 

Af^{t) = VAfo[ait) + a^t)] . (C4) 

Here [in the notation of Eqs. (2.34) - (2.37)] Wo = Aq^ is 
the carrier's photon flux and a{t), a^{t) are the time- 
domain side-band annihilation and creation operators 
with commutation relations [time-domain versions of 
(C2)] 

[a{t),a{t')]=Q, [a\t),aHt')]=0, 

[a{t),a^{t')] = 6{t-t') . (C5) 

It is straightforward, using these commutation relations, 
to verify that 

[M(t),M(0] = o. (C6) 

It is interesting to note that, although the photon num- 
ber flux self commutes, the energy flux (energy passing a 
fixed transverse surface per unit time) 

/■OO POO 

i{t) = n duj duj' \/^ ala^' (.i^ip-f ) ^^7) 
io io 



{£{t)Am^^- (C8) 

This can be thought of as due to the energy-time uncer- 
tainty relation for photons. On the other hand, when (as 
in gravitational-wave interferometers) the light consists 
of a monochromatic carrier plus signals encoded in side 
bands with frequency $7 = w — ^ oJq, then for all 
practical purposes, does self commute. 
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